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We investigate the impact of peculiar velocity effects due to the motion of the solar system 
relative to the microwave background (CMB) on high resolution CMB experiments. It is well 
known that on the largest angular scales the combined effects of Doppler shifts and aberration are 
important; the lowest Legendre multipoles of total intensity receive power from the large CMB 
monopole in transforming from the CMB frame. On small angular scales aberration dominates and 
is shown here to lead to significant distortions of the total intensity and polarization multipoles 
in transforming from the rest frame of the CMB to the frame of the solar system. We provide 
convenient analytic results for the distortions as series expansions in the relative velocity of the two 
frames, but at the highest resolutions a numerical quadrature is required. Although many of the 
high resolution multipoles themselves are severely distorted by the frame transformations, we show 
that their statistical properties distort by only an insignificant amount. Therefore, cosmological 
parameter estimation is insensitive to the transformation from the CMB frame (where theoretical 
predictions are calculated) to the rest frame of the experiment. 


I. INTRODUCTION 

The impressive advances being made in sensitivity and resolution of microwave background (CMB) experiments 
demand that careful attention be paid to potential systematic effects in the analysis pipeline. Such effects can arise 
from imperfect modelling of the instrument, e.g. approximations in modelling the beam a a or incomplete 

knowledge of the pointing, but also from more fundamental effects such as inaccurate separation of foregrounds (see 
e.g. Refs. ^ for reviews). In this paper we consider errors that may arise due to neglect of the peculiar motion of 
the experiment relative to the CMB rest frame (that frame in which the CMB dipole vanishes). For short duration 
experiments (e.g. balloon flights such as MAXIMA [7] and BOOMERANG [8]) the relative velocity is constant over 
the timescale of the experiment, but for experiments conducted over a few months or longer, and particularly for 
satellite surveys [9, 10, II], the variation in the relative velocity adds additional complications. In principle, the 
modulation of the aberration arising from any variation in the relative velocity must be accounted for with a more 
refined pointing model for the experiment [|l2| , [l3[ when making a map. 

For a relative speed of j3c (where c is the speed of light and (3 ^ 1.23 x 10“^ for the solar-system barycenter relative 
to the CMB frame), the r.m.s. photon Doppler shifts and deflection angles are /3/-\/3 and y^2/3/3 respectively. Despite 
these small values, significant distortions of the spherical multipoles of the total intensity and polarization helds do 
arise. A well known example is provided by the CMB dipole seen on earth, which, given the observed spectrum, 
arises from the transformation of the monopole in the CMB frame. More generally, on the largest angular scales 
the combined effects of Doppler shifts and aberration couple the total intensity monopole and dipole into the Zth 
multipoles at the level 0{f3^) and respectively. Given the size of the non-cosmological monopole, annual 

modulation of the dipole by the variation in the relative velocity of the earth in the CMB frame must be considered 
in long duration experiments. 

In this paper we concentrate on the effects of peculiar velocities on small angular scale features in the microwave 
sky. On such scales, aberration dominates the distortions and becomes particularly acute when the angular scales of 
interest, 0{l/l), drop below the r.m.s. deflection angle, i.e. I ^ 800 for the transformation from the CMB frame to 
that of the solar system. We provide simple analytic results for these distortions to the total intensity and polarization 
fields as power series in the relative velocity (i. The power series converge rather slowly at the highest multipoles for 
most values of the azimuthal index m [the leading-order corrections go like 0{l(3)] but the distortions can still easily be 
found semi-analytically with a one-dimensional quadrature. If the transformations of the multipoles carried through 
to their statistical properties, theoretical power spectra computed in linear theory (e.g. with standard Boltzmann 
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codes 011 ) would not accurately describe the statistics of the high resolution multipoles observed on earth. (The 
theoretical power spectra would still be accurate in the CMB frame.) It is straightforward to calculate the statistical 
correlations of the multipoles observed on earth. Fortunately, as we show here, the statistical corrections due to 
peculiar velocity effects turn out to be negligible despite the large corrections to the individual multipoles. It follows 
that for the purposes of high resolution power spectrum and parameter estimation, the transformation from the CMB 
frame can be neglected. 

This paper is arranged as follows. In Sec. II we describe the transformation laws for the total intensity multipoles 
in specific intensity and frequency-integrated forms. Convenient series expansions in (3 of the transformations are 
provided, and their properties under rotations of the reference frames are described. The statistical poperties of the 
transformed multipoles are investigated by constructing rotationally-invariant power spectrum estimators and full 
correlation matrices. In Sec. Ill we discuss the geometry of the frame transformations for linear polarization, and 


present power series expansions for the transformations of the multipoles. The behaviour under rotations and parity 
are also outlined. Power spectra estimators and correlation matrices are constructed, and cross correlations with the 
total intensity are considered. Some implications of our results for survey missions are discussed in Sec. IV, which is 
followed by our conclusions in Sec. 0. An appendix provides details of the evaluation of the multipole transformations 
as power series in (3. 

We use units with c = 1. 


II. TRANSFORMATION LAWS FOR TOTAL INTENSITY 


We consider the microwave sky as seen by two observers at the same event. Observer S is equipped with a comoving 
tetrad {(e^)“}, ^ = {0,1, 2,3}, and observer S' carries the Lorentz-boosted tetrad {(e(^)“}. The relative velocity of S' 
as seen by S has components on {(ei)“}, i = {1, 2,3}, which we denote by the spatial vector v, which has magnitude 
(3. The S observer receives a photon with four-momentum p°‘ when their line of sight is along n, so the photon 
propagation direction is —n. For S the photon frequency is v where hv = Pa{eoY Planck’s constant), while S' 
observes frequency 


v' = - 1 - ri • v). 


where 7 ^ = I — . The line of sight in S' is 


"h-v-|-/3\ .. n—n - vv 
n' = I -: —- 1 V -I- 


1 -I- n • V 


7(1 -I- n • v) ’ 


( 1 ) 


( 2 ) 


where v is a unit vector in the direction of the relative velocity. 

Denoting the sky brightness in total intensity seen by S as /(j^, h), the brightness seen by S' is (e.g. Ref. |M) 


/'(i/',h') = /(zz,h) 


(3) 


If S and S' use their spatial triads {( 6 ^)“} and {(e')“} to define polar coordinates in the usual manner, and expand 
the sky brightness in terms of scalar spherical harmonics, i.e. I{v,n) = we find the following 

transformation law for the brightness multipoles: 


= 0 / dn7(l + n-v)a/,^,(j/)Wm'(n)Ri™(n'), 


(4) 


where v = + n • v)“^, and we have used i/'^dri' = i/^dn. 

It will prove more convenient to consider the integral of the brightness over frequency, / (n) = J dv I{v, n) . The 
transformation law for this flux per solid angle follows from integrating Eq. (^: 

(5) 


/'(n') = /(n), ^ 

Expanding /(h) in spherical harmonics, we find the multipole transformation law 

aL = 04m' [ dh[7(l-hh- v)]2Wm'(h)Y,;,(h'), 
— 'y ' Ar(im)(im)'Q;'m' 


( 6 ) 
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where = f dvaj^^v). The second equality defines the kernel which relates the frequency-integrated 

multipoles in S and S'. Dividing by four times the average flux per solid angle gives the multipoles of the 
gauge-invariant temperature anisotropy in linear theory (e.g. Refs. 0 @)- 

If we choose the spacelike vectors of the tetrad {(e^)“} so that the relative velocity is along ( 03 )“, the multipole 
transformation law becomes block-diagonal, oc 5mm', with no coupling between different m modes. The 

kern el fo r a general configuration can then be inferred from its transformation properties under rotations described in 
Sec. |II A|. In the appendix we evaluate Eq. (0) as a series expansion in (3 for general spin-weight functions, including 


terms up to 0(/3^), for the case where v is aligned with ( 63 )“. The expression is cumbersome, partly due to the fact 
that the transformation law is non-local in frequency. For large I the aberration effect dominates Doppler shifts and 
the frequency spectrum of the multipoles is pr eser ved by the transformation. We also give the result obtained by 
integrating over frequency; setting s = 0 in Eq. (^) we find the series expansion of the kernel K^im){imY up to 0(/3^): 


K(im)(i'm) = 5w ^ — 1){1 — 2) + Cfm{l + 2){l + 3) + — l{l + 1) + 2^ 

+ 5i(^ii + i)PClmil + 3) — — 2) 

+ 2 -I- 3) -I- dz(i'-2)/3^C'(Z-|-2)mC'(i+l)m 2 ~ ~ 2), 


where Cim = oC'/m with 




Im — 


'(/2-m2)(F_s2) 

12 ( 4/2 _ 1 ) 


(7) 


( 8 ) 


Comparison with Eq. (A 6 ) shows that for high I the aberration effect described by the term in Eqs. (Q) 

and is dominant. For /31 > 1 the series is slow to converge for \m\ <C I since the leading-order corrections go 
like 0{ip), reflecting the fact that the deflection angle due to aberration is comparable to the angular scale of the 
spherical harmonics at this 1. For /3 « 1.23 x 10“^, appropriate for the solar-system barycenter relative to the CMB 
frame, /31 ~ 1 corresponds to multipoles I ~ 800. In this case, the kernel K^im)(im)' is easily evaluated by a numerical 
quadrature. We show some representative elements of the kernel in Fig. |^, which demonstrates that the multipoles do 
indeed suffer severe distortion for 1 > 1//3, as suggested by Eq. (^. For \m\ close to I, sCim ~ 0(1“^/^) and retaining 
only the terms given in Eq. § is accurate to much better than 0.1 per cent for the I range probed by e.g. Planck 
(1 ^ 2000). For such values of m the distortions to the multipoles are only small, with leading-order corrections at 
1' = 1 ± 1 of 0{P'/l)- For ip 1, the departures of the kernel from the identity SwSmm' are very small, giving 
negligible distortions to the multipoles except for I close to unity when the non-zero coupling to the (large) monopole 
can give significant distortions, as described in Sec. I. 


A. Rotational properties 

If we rotate the relative velocity v to keeping the tetrad (e^)“ fixed, [thus inducing a transformation of the 
Lorentz-boosted tetrad (ej^)®], the frequency-integrated multipoles continue to be given by Eq. (^, but with v replaced 
by Z?v in n' [Eq. (^] and in (1 -|- n • v). With the change of integration variable n ^ Dn, the integral defining the 
transformed kernel Rl(im)(im)'(TIv) becomes 

J dh"f'^{l + n-w)Yiim'iDn)Yi*^{Dh') = J dh^'^{l + h- v)D~^Yi,m'iii-)[D~^Yim{n')]*, (9) 

where D~^Yiim'{n) = Sm" ^m'm"^irn" {n) with a Wigner D-matrix. (Our conventions for D-matrices follow 

Refs. [^ ^.) It follows that the transformed kernel is given by 

K(lm)ilmy{D'v) = ^ (10) 

MM' 

Instead of rotating the (physical) relative velocity of S and S', we could imagine rotating the spatial triad (e^)® ^ 
D{ei)°‘. Under this coordinate transformation, the Lorentz-boosted frame vectors transform similarly: (e')“ ^ Z?(e')“. 


^ Here, D denotes the appropriate representation of the rotation group 5'0(3). 
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FIG. 1: Representative elements of the frequency-integrated kernel evaluated with the relative velocity (/3 = 1.23 x 

10“^) along|_(e 3 )“. The results of a numerical integration of Eq. are shown in dark gray, while results based on the series 
expansion ([^ are shown in light gray. The smaller (absolute values) of the two are shown in the foreground. Elements are 
shown for I = 1500 (left), I = 700 (middle), and I = 50 (right), with m = 0 (top) and m = I (bottom). 


For a fixed sky, the multipoles seen by S and S' transform according to e.g. ^ Sm' (which is equivalent 

to rotating the sky with D~^ leaving the tetrad fixed). It follows that under coordinate rotations, the kernel transforms 
as 






^^M'rn' 


( 11 ) 


Note that the (passive) rotation of the frame vectors by D~^ has the same effect on the kernel as the (active) rotation 
of the relative velocity v by U, as expected. 

Finally, we consider (active) parity transformations v ^ — v with the tetrad (e^)“ held fixed. Using ¥)„(—n) = 
i-iyYimin) it is straightforward to show that 




K^lrnWmyM- 


( 12 ) 


The behaviour of the kernel under parity ensures that if we simultaneously invert v and the sky [af^ ^ (— l)^af^], 
the multipoles seen by S' transform to (—l)^afm. 

These transformation properties of the kernel under rotations allow one to generalise Eq. (^ easily to the case 
where v is not aligned with ( 63 )“. 


B. Power spectrum estimators 

We have seen how aberration effects lead to significant distortions of some of the high-/ multipoles in transforming 
from the CMB frame to the frame of the experiment. In the next two subsections we investigate the impact of these 
distortions on the statistical properties of the multipoles. 

We assume that in the CMB frame (S') the second-order statistics of the anisotropies are summarised by 


= Cl^SwSmm', 


( 13 ) 
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FIG. 2: Representative elements of the kernel Wu/ evaluated with relative velocity (3 = 1.23 x 10 The results of a numerical 
integration of are shown in dark gray, while results based on the series expansion ( 0 ) are shown in light gray. The smaller 
(absolute values) of the two are shown in the foreground. Elements are shown for I — 1500 (left), I = 700 (middle), and I = 50 
(right). 


appropriate to a statistically-isotropic ensemble with power spectrum . (The averaging is over an ensemble of 
CMB realizations.) It is this C/ for ^ > 2 that is computed with linear perturbation theory in standard Boltzmann 
codes (e.g. Refs. 0E1)- 

We begin by considering the quadratic statistic 

(14) 

m 

which is evaluated by S'. In the absence of noise this statistic is the optimal (minimum-variance) estimator for the 
power spectrum if we ignore peculiar velocity effects. By construction, C^' is independent of the choice of spatial 
triad, but is only invariant under rotations of the sky in S (af^ ^ relative velocity is also rotated to 

Hv. However, averaging over CMB realizations keeping the relative velocity fixed we obtain a quantity {Cl^') which 
is obviously invariant under rotations of the sky in S. The average which determines the bias of the power 

spectrum estimator C^^', is linearly related to the 

I'mm' 

= Y^Wu'C''. (15) 

The kernel Ww depends only on the relative speed (3 and not the direction go we can always evaluate it with v 
aligned with (eg)®. 

The series expansion (|^) of can be used to evaluate Ww. Correct to we find 

Ww = Sw (l - + 1-8)^+ 3{2l + l)^^- 

Again the series is slow to converge for ^ > 1 and the terms neglected in Eq. (0) are non-negligible. We show Ww 
for some representative I values in Fig. g. It is clear from the figure that Ww is well localised in comparison to any 
features in the CMB power spectrum for the range of (3 of interest here. In this case, we can approximate 

(A"') « 

I' 

= C/^[1 + 4/32 + 0(/33)]. (17) 

The effect of the velocity transformation is thus to rescale the amplitude of the power spectrum by 1 -1-4/3^. This bias 
is clearly insignificant. is actually independent of I to all orders in (3. To see this we form Wi' directly 
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using the integral expression (H) for K(irn){im)'- The result simplifies to 

= J dh (1 + ii • v)^ = 7^ 2/3^ + i/3^^ , (18) 

on using the completeness relation 

sYim{ni)sYi*^{n2) = (5(ni - fiz), (19) 

Im 

and the addition theorem for the spherical harmonics. The series expansion of Eq. ( |T8| ) agrees with Eq. (|l7|). We 
conclude that despite the fact that the multipoles themselves can be severely distorted by aberration for 1 in 

passing from the CMB frame to that of the solar system, the quadratic power spectrum estimator is negligibly biased 
since the effect of the velocity transformation is to convolve the power spectrum with a narrow kernel Wui that sums 
to very nearly unity. 


C. Signal covariance matrix 


Assuming Gaussian statistics in the CMB frame, the multipoles in S' will also be distributed according to a 
multivariate Gaussian since the transformation (i is linear. In this case, the covariance matrix (af^af/m/) contains 
all statistical information about the anisotropies in S', and as such is an essential element of optimal power spectrum 
estimation. 

If we make use of Eq. (H) , the covariance matrix in S' reduces to 




LM 






( 20 ) 


The presence of the preferred direction v breaks statistical isotropy in S', and the multipoles are correlated for I ^ I' 
and m ^ m'. The structure of the covariance matrix in S' depends on the choice of the spatial triad (e^)® with respect 
to the relative velocity of the two observers. Aligning ( 63 )“ with v, the m-modes decouple in and so also 

in the covariance matrix. Furthermore, for the values of /3 of interest here (/3 <SC I), the kernel falls rapidly 

to zero for I and I' differing by more than a few (see Fig. |^), so the same will be true of the covariance matrix. It 
follows that we can approximate 


(afmaf'm' 


LM 


{lm)(LM) Ylm,)' {LM)' 


( 21 ) 


(Pulling out C[/ instead will give essentially the same result for smooth power spectra.) The summation in Eq. ( [^ ) 
is most easily evaluated by substituting the integral representation (||) for Arpm)(LM) and using the completeness 
relation (p^. We find that 'Llm reduces to the integral 

J dn [7(1 + n • v)]*Yi*^{h')Yi,m'{n') = J dn' [7(1 - n' • v)]-^Yi*^{h')Yi,m'{n'), (22) 

where we changed the integration variable to h' and used 7(1 + n • v) = [7(1 — h' • v)]“^. Note that both spherical 
harmonics have the same argument in the integrand, so we don’t expect the same 0{ip) terms at high I that arise in 
the kernel Kty^y^yi^y^y. Eq. ( p 2 | ) can easily be evaluated for ( 63 )“ along v (in which case there is no coupling between 
different m) by expanding in /3: 

E “ 1)] + ^Z(i' + l)6/3C'im + 5z(i'_l)6/3C'(;+i)m 

LM 

+ <5i(i'+2)21/3^GimC'(;_i)m + <5;(i/_2)21/3^C'(i+2)mC'(i+l)m + O{0^). (23) 

This result for the covariance matrix in S' could easily be used in maximum-likelihood power spectrum estimation 
(see e.g. Ref. [^) to correct for the bias due to peculiar velocity effects. However, since the leading corrections are 
only 0{/3), even at high I, the effects will be negligible. 
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III. TRANSFORMATION LAWS FOR LINEAR POLARIZATION 


The linearly polarized brightness in S is described by Stokes parameters Q{v, n) and U{v, n). The Stokes parameters 
depend on a specific choice of orthonormal basis vectors {mi, m 2 } for each line of sight n. If {mi,m 2 ,— n} form a 
right-handed orthonormal set, the Stokes parameters are related to the linear polarization tensor by 


V 


ab 


-[Q(mi (g) mi 


m2 (g m2) -I- C/(mi (g m2 -I- m2 (g mi)]. 


(24) 


The Stokes parameters transform under changes of frame in the same way as the total intensity, i.e. 

/ ^ \ ^ 

Q'{iy\h') = Q{iy,h) , 


(25) 


and similarly for t/, provided that the basis vectors are transformed according to 

m' = mi -I- (7 - l)mi • vv - ym^ • vn'. 


(26) 


where i = 1,2. It is straightforward to verify that this transformation law preserves orthonormality, and also that m' 
is obtained from m^ by parallel transport on the unit sphere along the great circle through n and n' (and so through 
V also). In terms of the polarization tensor, the frame transformation law can be written as 


iP'“'’(:/',n') =iP|f(i/,n;v) 



(27) 


where n; v) is h) parallel propagated to n'. The 1-1-3 covariant form of this transformation was given in 

Ref. ||. 

If S and S' introduce polar coordinates as in Sec. II, the polarization tensor can be expanded in symmetric trace-free 
tensor harmonics 

Vab{v,il) = i X! («): (28) 

Im 

which dehnes the electric (E) and magnetic {B) multipoles. Using Eq. (^^ we can extract the multipoles seen by S'. 
For we find 


E/ 

^Im 


{E) = 


dn7(l 


n • V a 


)[° 






(29) 


with a similar result for Here v) is parallel propagated to n', and similarly for the curl 

harmonics. Note how in general the frame transformation mixes E and B polarization. Equation (|^) is valid quite 
generally, and is useful for discussing the rotational properties of the transformations (see later). However, to compute 
the transformation laws it is again convenient to arrange (ei)“ so that v is along ( 63 )“. We can then exploit the fact 
that the polar basis vector helds 0 (n) and (f>{h) are parallel propagated along longitudes to simplify (n; v). The 

gradient and curl harmonics can be written in terms of spin-weight ±2 harmonics (our conventions follow Refs. ||^, p^): 


\^G ab 
^ Im 


ab 

^ Im 


—^{-2Yimriv (g m -I- 2YimrrY (g m*), 

v 2 

—j={-2Yirara (g m - 2Eimm* (g m*), 

iv2 


(30) 

(31) 


where the complex vector m = (0 -|- z0)/-\/2, so that Eq. (j^) can be written as 


{af:,±^aZ){E 


) = ±iaf^,){y) J dhj{l + n ■ v)± 2 Yi,^>{n)± 2 Yi*^{h'). 


(32) 
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Our and are \/2 times the gradient (G) and curl (C) multipoles introduced in Refs. | |23[ With this convention the power 
spectra of the electric and magnetic multipoles agree with those defined in the spin-weight formalism [Ea, Gfil . 



The integral on the right-hand side is evaluated as a power series in (3 for general spin-weight s in the appendix. 
For our purposes it will be more convenient to consider the frequency-integrated multipoles, e.g. af^ = 


Integrating Eq. (29) over frequency, we find 


I'm' 

B/ ^ Jy' B • jy E 

^Im / V -\--^ilm){lm)'^l'm' — ■> 

I'm' 


(33) 

(34) 


where the kernels 


-\-^{lm){lm)' 
— ^{lm){lm)' 


I dn[7(l-hn.v)]2l^«)),V(“;'^)^(t)a.(A')= / 

-i J dn[ 7 (l-f n-v)] 2 l^|'(,^^),(n;v)Y(f^)„(,(n') = 


dn [ 7(1 + n • (fi; v)Y(f^)„j,(ri'), (35) 

i J dri[ 7 (l-f n-v)]2Y|G^^),(h;v)y(f^)^(,(n').(36) 


The behaviour of ±K(^im)(im)' under rotations v ^ Dv is the same as for the total intensity kernel, Eq. (H), since 
the tensor harmonics transform under rigid rotations with the same D-matrices as the scalar harmonics ^ . This 
property of the tensor harmonics also ensures that under rotations of the coordinate system, (ei)“ ^ Z?(ei)“, the 
electric and magnetic multipoles transform irreducibly to e.g. Em' ^^m'm'^fm'- Under inversion of v with (e^)“ held 
fixed, the kernels transform to 


+ ^{lm){lm)' 


(-V) = (-1) 


l+l' 


+ ^{lm){lm)' 


(v), 


-^{lm)(lm)'{ v) — ( 1 ) —-^(/m)(/m)'(^); 


(37) 

(38) 


so that under simultaneous inversion of the sky in S', af^ (-i)'«L and af^ inversion of v, the 

multipoles in S' transform like those in S. 

The frequency-integrated kernels are most simply evaluated with v alon g (e .^)”. In this case the m-modes decouple, 
as with the total intensity. Writing ±K = { 2 K ± _2if)/2, we can use Eq. {M), which evaluates sKiy^nyi^y as a series 
correct to 0{(3'^), to show that 


+ ^(lm)(l'm) 


= Sw 


1 


— I)(Z — 2) -I- 2Qm(^ + 2)(^ -I- 3) -I- TO — ;(Z -I-1) -I- 6 — 


4to^ 


-b 


24to^ 


i{i + i) p{i + iy 


+ Sl{l' + l)P2Clmil + 3 ) — ~ 2 ) 

+ ^l{l'+2)P^2Clm2C(^i_i')jn-{l + 2){l -|- 3) -|- -2)/3^2C'(/+2)m2C'(i+i)m - (Z — 1)(Z — 2), 


(39) 


and 


-K, 


&f3m 


6/3^to 


6/3^to 


{lm)(l'm) — 1(^1 _y <5j(Z' + l) 2 C'im(Z -b 3) ^ ^ -b l5i(i/_i)2C'(;+i)m(Z 2 ) ' 


(40) 


(Equivalent results, correct to 0(/3), have already been worked out in l-b3 covariant form ||2^.) The kernel -KtimVimV 
is suppressed at high 1. It receives comparable contributions from Doppler and aberration effects for all I [see Eq. (|A(j )] 
in contrast to the +K(^i^yi^y and the total intensity kernel 7^(im)(im)' which are dominated by aberration effects at 
high 1. The series expansion of +K(^im){imY is slow to converge for Z/3 ^ 1 when |to| ^ I, and there are large distortions 
to the electric and magnetic multipoles for these indices. Electric multipoles nearby in I couple in strongly to distort 
aJE} S'lid similarly for the magnetic multipoles. Eor I ^ 1 the kernel +K(^im){imy is almost indistinguishable from the 
total intensity kernel Kyi^^yi^^y. The cross contamination of e.g. B hy E due to the frame transformation is much 
weaker, with the maximal effect ^ 0{(3/l) at leading order occuring for |to| « 1. [Note that, as with +K(^i^yij^y, the 
convergence of of Eq. ( p0[ ) is slow for Z/3 > 1 when |to| ^ L] The transfer of power from E to B is potentially the 
most interesting effect since in the absence of astrophysical foregrounds, inflationary models predict that magnetic 
polarization in the CMB frame on scales larger than a degree or so arises only from gravitational waves. However, on 
these scales a gravity wave background comprising only one percent of the large-angle temperature anisotropy would 
have B power far in excess of that generated in the frame of the experiment by transforming E from the CMB frame. 
On sub-degrees scales, where any primordial i?-polarization is expected to be very small, other non-linear effects, 
most notably weak lensing of E |28|[, will dominate the B signal produced by the velocity transformation. 
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A. Power spectrum estimators 

The second-order statistics of the polarization multipoles in the CMB frame, assuming statistical isotropy and 
parity invariance, define power spectra: 



= SlliSmm'Cf'^, 

(41) 


= SwSmm'CP'^, 

(42) 

Slra^l*m') 

= SwSmm'Ci^, 

(43) 


with no correlations between B and E or I. We can form estimators of these power spectra from the multipoles in S' 
by analogy with Eq. (0), e.g. 

= (44) 

m 

Since these estimators are rotationally invariant we can compute them for v aligned with ( 63 )“ using Eqs. and 
( ^ ) without loss of generality. 

The expected values of the power spectra estimators can be expressed in terms of the power spectra in the CMB 
frame using Eqs. (^), (p^, and (p^): 


//^EEf\ 

4 

\+K(lm)ilmy\'^Cl?^ + \-K(im){lmy\^Cl^^, 

(45) 

\^l / 

21-1-1 f—f 

1 m m 

{CD 

4 

\+K(lm)ilmy\‘^CP^ + \-K(im)(lmy\^C^^, 

(46) 

21 + 1 ^ 
I'm'm 

{CD 

4 y^ 


(47) 

21-1-1 

1 m m 


Substituting the power series expressions for the kernels and performing the summations over m 


and m' we find 


and 


1 

21 + 1 


^ ( \ + K{lm){lm)' I 

mm' 


1 

2r+T 


E 

mm' 


-K, 


{lm){lm)' \ 


(1 - + 4 )(i - 3 )) + 




(l-2)^(/ + 3)(Z- 1) 
3(/ + 1) (2/ + 1) 




12 


1 

21 + 1 


^ ^ +^{lm)(lm)'^(lm){lm)' 
mm' 


&. (1 - + / - 10)) + 

+ 4«,.-.)/3V(' + 3)('- 1 ) 3 ) 2 ;+’l)' 


(48) 

(49) 


(50) 


correct to 0{f3'^). Eor I 1 the right-hand sides of Eqs. (0 ) and (0 are almost equal to each other and to the kernel 
Will which determines the bias in the total-intensity estimator C[' . As with the total intensity, the series in Eqs. 0 - 
p0|) are slow to converge for 1(3 > 1. The bias of by A-polarization is controlled by J2mm’ P/(2^ +1), 

which falls off rapidly with 1. In Fig. || we compare this contribution to the expected with the B-polarization 

power spectrum due to primordial gravity waves and weak leasing of the E-polarization. The cosmological model is 
a Lambda, cold dark matter (ACDM) in which gravity waves contribute one percent to the large-angle temperature 
anisotropy. As remarked earlier, the contamination arising from the frame transformation is well below the expected 
in such a model. 

The means of the estimators and , defined by analogy with would vanish in the absence of peculiar 

velocity effects (and foregrounds) due to parity. The velocity transformations preserve these zero means since 


mm' mm' 


(51) 
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FIG. 3: Contribution of to the mean estimator {Cp^') in a ACDM model with one percent contribution to the total- 
intensity quadrupole from gravity waves. This velocity effect is compared with (in the CMB frame) due to primordial 

gravity waves (solid line) and weak lensing of the if-polarization (dashed line). 


These results are easily proved by choosing v along ( 63 )“ so that all kernels are real, and using the general results 

^{l-Tn){l' -m')- 

fall off sufficiently rapidly with — ^| t hat t hey are 




^ {lrri){lrri)' 

The kernels represented by the left-hand sides of Eqs. 


we can 


narrow compared to expected features in the primordial power spectra Following the analysis in Sec. IIB 
pull out Cy aX V = I from the summations in Eqs. (E^-(^). Performing the sums over V, we find 


1 


21 + 1 


1 


21 + 1 


E 1+^' 

mm'l' 




= 1 + 4/3- 


l^ + l-3 

l{l + l) ’ 




= 


12 


mm'l' 


i{i + i)’ 


1 


21 + 1 




mm'l' 
2\ 






(52) 

(53) 
, (54) 


correct to 0(/3^). For large I the right-hand sides of Eqs. (|^ ) and ( ^ approach 1 + 4/3^; as with the total intensity, 
there is a negligible scaling of the amplitude of the power spectra estimated in the S' frame due to the frame 
transformation. Note that 


1 


2 / + 1 


E ii+A'i 


(lm){lm)' \ 


-K, 


1 


{lm){lm)' \ 


mm'l' 


2 ( 2 / + 1 ) 


E 


{lm){lm)' 




mm'l' 


= 7" (^1 + 2/3^ + -/3^ 

where we have used the completeness relation and addition theorem for the spin-s harmonics. Adding Eqs. 
(m) we obtain the series expansion of the exact result in Eq. (M). 


(55) 

and 


^ At low I the polarization power spectra vary rapidly (as power laws) with 1. Over this part of the spectrum the approximation that the 
power is approximately constant over the width of the convolving kernel is still valid since the latter are essentially Kronecker deltas at 
low 1. 
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B. Signal covariance matrices 


The calcul ation of the covariance matrix of the polarization multipoles in S' follows that for the total intensity 
given in Sec. IIC. For smooth power spectra we can approximate 


{aZ4l') 


E + + cr E 

LM LM 

cr E +^(^-)(LM) + E^m)'(LM) + cr E 

LM LM 

E +^(l'm)(LM)K*l'm)' {LM) ■ 


{lm){LM)-K*lmY {LM)-! 

(56) 

-~7<^{tm){LM)-K*lm)'{LM)^ 

(57) 


(58) 


LM 


The remaining correlators would vanish for v = 0 due to parity invariance. For non-zero v we can approximate 


aZaf^Z) ^ 

S iCM E +^(lm)(LM)-Kamy(LM) 



LM 

LM 

aZ4l') ^ 

S —iC/^ E -^(lm)(LM)K*im)'(LM)- 



LM 


(59) 

(60) 


If we align v with ( 63 )“ we can evaluate these expressions by substituting for the series expansions of the kernels from 
Eqs. (1^) and (|4^). The m modes decouple and we find 


LM 


LM 


1 + 3P^ 


LCf,+l)r. 


72 C, 


bn 


16m^ 


-1 


+ 5i(ii + i)%f32Clm + l^;(i'-l)6/32C'(Z-|-l)m 
+ <5i(i'+2)21/3^2C'im2C'(i-l)m + <5i(i'-2)21/3^2C'(;+2)m2C'(Z-|-l)r 
36/3^m^ 


(61) 

(62) 


E +^(lm){LM)K{l'm){LM) — + 1)(^ + + ^ll+l)m) ~ ~ 18)(2C'fm + Qm) 

LM ^ 

-1-18-1- 2(1 — 2)^2C(/+l)mC'(/+i)m + 2(^ -f 3)^2C'imC/m]^ 

~ 5l{l' + l)P[(^ — 3)Clm — (^ + 3)2C'/m] + 5i{ii-i)f3[(l + 4)C(i+i)m — (I — 2)2C'(/_|_i)j„] 

+ 2^l{l'+2)/3‘^[il + 2)(/ -I- 3)2Cim2C{l-l)m + (^ ~ 3)(1 — 4)C'imC'(Z-l)m 

— 2(1 + 3)(l — ‘^)2C(l-l)mClm\ + 2^l{l'-2)P'‘^[(^ ~ 1)(^ ~ 2)2C'(;+2)m2C'(i-|-l)m 
+ (^ + 4)(^ + 5)C(/_|_2)mC'(Z-|-l)m ~ 2(Z — 2)(l + 5)2C'(i-|-l)mC'(i-|-2)m]i (63) 


correct to 0(pS). This final expression is cumbersome and hides the fact that the leading order corrections to the 
covariance matrices are only 0(/3), rather than 0((3l). To see this, we can expand Eq. (|^) in 1/1 for large I to find 


E +^{lm){LM)K*l'm){LM) 
LM 


Sw 


1 + /3 


15 155 -k 84m2 


8P 


1 3(7 + 4m2) 

+ V+i)/3(3-y - 


+ ^i{i'-i)P ( 3 + y ~ 


1 29-kl2m^\ , , ^2 /"21 5 159-k 84m" 


8F 


+ v + 2)/3 [^-jr 


16/2 


+ V-2)/3 1^ + ^- 


21 7 223 -k 84m" 


16/2 


(64) 


correct to 0(1 2). For |m| « / the expansion in 1/1 is slow to converge, and the full expression, Eq. (|^, should be 
evaluated exactly if the (very small) corrections to the covariance matrices are to be included in a statistical analysis. 
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It is worth noting that 




LM 


j dn[ 7 (l-n-v)] ®[ 2 y;;,(n) 2 r/',„(n) + _ 2 yi;;,(n)_ 2 yi'™(n)], 


(65) 


for V along ( 63 )“, where we have used the completeness relation, Eq. (|19[). It is straightforward to show with an 
expansion in f3 that Eq. (p^is consistent with adding Eqs. and (p^). 

Eor the correlators and (aj^aj/^,), which would vanish for v = 0, we require the results (for v aligned 

with ( 63 )“) 




(lm)(LM)-K*l'm){LM) 




LM 


l{l + l) 




— ^i(r-i)2C'(i+i)m(7^ + 4)- 






(lm){LM)K’^l' ra){LM) 


LM 


l{l + l){l + 2y 

- <' - IX' - + 

6 0 ‘^Tn 

+ + l)(^ - 2)2C'(;+i)m - {I + 4)(l + 2)C{i+i)m\ _|_ 2 ) > 


correct to 0{P‘^), and the general result 


E+ E +^(l' -m'){LM)-K*i -Tn)(LM)- 


( 66 ) 


(67) 


( 68 ) 


LM 


LM 


The leading order corrections to the components of the correlation matrices that vanish for v = 0 are 0{(3), and are 
suppressed at large 1. and so can safely be ignored. Eor completeness we note that 






LM 


\ J dh[ 7 (l-n-v)] ®[ 2 yi;;,(n) 2 Ep™(n) - _ 2 l));,(n)_ 2 yi',„(n)]. 


(69) 


This result is easily shown to be consistent with Eqs. (@ and (p§). 


IV. IMPLICATIONS FOR SURVEY MISSIONS 

For experiments which observe for less than a month or so the velocity of the instrument relative to the CMB 
frame can reasonably be considered constant. In this case a map in the frame of the instrument can be made with 
no account of the effects considered in this paper. Accounting for the peculiar velocity relative to the CMB frame 
can be deferred until the statistical properties of the map are considered. As we have shown here, peculiar velocity 
effects can safely be ignored when estimating smooth power spectra since the estimated power spectra are essentially 
convolutions of the spectra in the CMB frame (which we can reliably compute with linear perturbation theory) with 
narrow kernels that integrates to unity. 

For survey experiments that observe for the order of a year or more the variation in the orbital velocity of the 
instrument adds another potential complication. Modulation of the dipole by the orbital velocity of the earth was 
visible in the COBE DMR data here we are interested in effects at small angular scales. To estimate the 
importance of the effect we consider a toy model of the Planck High Frequency Instrument (HFI). We approximate 
the orbit of the satellite relative to the sun as a linear motion with [3 = 10“'^ for six months, after which the direction 
of motion is reversed for the next six months of observation. Clearly, this toy model will over estimate the effects of 
the variation in orbital velocity. Planck will cover the full sky in six months, so for each six month period we could 
make a map and extract the spherical multipoles. In our toy model these two maps are produced in frames with 
a relative velocity of 2/3 = 2 x 10“^. In the 1-range relevant to Planck we need only retain the 0(/3) corrections in 
Eq. ( 0 ), so the difference between the multipoles measured from the two maps can be approximated as 




( 70 ) 
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for large 1. Here, are the total intensity multipoles in the rest frame of the solar system. The r.m.s. difference in 
the multipoles is 


(|AaL|2)V2 « V2l3Wl-myP^ < V2f3l^, (71) 

which should be compared to the instrument noise. For the 100 GHz Planck HFI channel, the one-year pixel noise 
is 6.0 /iK in 9.2 arcmin (the beam full-width at half maximum) pixels. The noise on our six month maps will be 
larger than this figure by -\/2. A comparison of the noise on the recovered multipoles with the r.m.s. error due to 
the difference in orbital velocity shows that the latter is just above the noise in the region of the first acoustic peak 
in Cl^ (at I ^ 200) for \m\ small compared to 1. Combining maps at different frequency would reduce the noise 
while preserving the peculiar velocity effect. However, since we have certainly over estimated the importance of the 
variation in orbital velocity, it is likely that the variation in aberration due to the orbital motion of the earth need 
not be considered beyond the dipole (which is modulated by the large CMB monopole). In principle, the modulation 
of the high I multipoles could easily be accounted for during map-making by including the aberration corrections in 
the pointing model of the instrument . 


V. CONCLUSION 

We have shown that for total intensity the effect of the frame transformation from the CMB frame to that of the 
solar system produces large distortions in certain multipoles at high 1. These effects arise principally from aberration 
rather than Doppler shifts. The linear polarization multipoles are similarly distorted at high Z, but with the additional 
complication that there is some transfer of power between E and B polarization. This transfer is suppressed at large 
I, and receives comparable contributions from aberration and Doppler shifts on all scales. Although the power in 
B polarization is expected to be much smaller than that in E in the absence of foregrounds, the B polarization 
generated from E is well below the primordial level even if gravity waves contribute only one percent of the large- 
angle temperature anisotropies. If the gravity wave background is much below this level, weak gravitational lensing 
will dominate the primordial signal on all scales. This lensing signal is expected to be an order of magnitude larger 
than the B polarization generated from the frame transformation on large scales. 

Despite significant 0{f3l) distortions of certain multipoles at large I, peculiar velocity effects are suppressed in power 
spectrum estimators and the covariance matrices for the CMB signals. The effect of the frame transformation on 
the mean of the simplest power spectrum estimator is to convolve the spectrum in the CMB frame (which we can 
compute reliably with linear perturbation theory) with a narrow kernel that integrates to unity. For smooth spectra 
there is negligible bias introduced by such a convolution. For linear polarization, the bias of e.g. the H-polarization 
power spectrum by E is suppressed at large I, and is expected to be negligible on all scales. We also showed that 
the frame transformation has only a negligible effect [0(/3) as opposed to 0{f3l)] on the signal covariance matrices for 
smooth underlying power spectra. The leading order effect is a coupling to the adjacent I values, Z ± 1. For linear 
polarization additional correlations are induced between B and E polarization, and B and total intensity 7, since the 
frame transformation does not preserve parity invariance, but their level is negligible. 

If the CMB fluctuations are Gaussian in the CMB frame, the multipoles will remain Gaussian distributed in any 
other frame since the transformation is linear in the signal. The transformation does break rotational and parity 
invariance, however, and so the aberration effects described here may be important when searching for weak lensing 
effects in the microwave background (using small patches of the sky over a coherence area of the weak shear), or the 
effects of non-trivial topologies. 
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APPENDIX A: SERIES EXPANSION OF THE MULTIPOLE TRANSFORMATION LAWS 

In this appendix we outline the evaluation of the transformation law for the brightness multipoles as a power series 
in /3. We align the relative velocity with the vector ( 63 )“ so that there is no coupling between different m modes. To 
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allow us to discuss both total intensity and linear polarization, we consider the integral 


aimW)=Y^ f dh—ai'^{iy),Yi,^{n)sYi*^{n'), 
V '' ^ 


/2 


d 2 


where v'/v = 7(1 + n • v) and h' is given by Eq. (||). We Taylor expand ai'm{v) as 

ai'm{v) = + 

where /i = n • v, and we handle sYimi^') with the expansion 


2 dz /'2 


avm{v')r^ + 0(/3 ), 


= sYvm{d) - /3(1 - - 1 ) + o{P^). 

The derivatives with respect to /i can be eliminated with repeated use of the identity [po[ 

(/^ ~ ^sC(i+l)ms^(i+l)m T 2^^ sYlm T ^)sClmsY(^i_i'j^, 

where sCim is defined in Eq. (||), and residual factors of /r can be absorbed with the identity 

sm 

IJ'sYlm — sC{l+l)msY(^ij.i'jm 1(^1 sYlm T sClmsY(l—l)rn- 
With these results, we find the following expression for 


(Al) 


(A2) 


(A3) 


(A4) 


(AS) 


= i 1 - 


Psm 

l{l + l) 


2 - 1 /' 


/ d 
du' 




+ + ^)(^ + 2 ) - 2 (/ + + 2 + 1 ) + 1 - 

, d d^ y 


2 2 
s m 


2 2 

s m 


2l{l + l) Pil + l)"^ 


1 - j/' 


di/' 


' di/'2 y 






— fisCln 


-{I + 2 ) + u'— - 


dv' l{l + 2) 
Psm 


2(Z - 1) - - 2)i/'-y. _ 


d 

diy'^ 
d 


-2{l + 2) + il + 3y^-u'^ — 




iy ) 


+ -^0^sClmsC{l-i)m (^{1 + 1)(^ + 2) — 2,(1 + 1)^^ ^ 0‘[l-2)my) + 0(y ). (A 6 ) 

Integrating this result with respect to v', the kernel sK(^im)(imy introduced in Sec. II (s = 0) and Sec. (s = ±2) 
evaluates to 




1 — T7^ -TT + -y (sCn.^^yi — !)(/ — 2 ) + sCfyi + 2){l + 3) 


l{l + l) 2' 


+ ly + — i{i +1) + 2 — 


2 2 
s m 


6s^m^ 
1(1 + 1 ) ^ 12(1 + 1)2 
3/3sm 


+ + 3)(^1- -2)(^1- 


: + i)(i- 

l-2)(l + 3 

with K(im){imY = 0 for m ^ m' in the configuration with v along ( 63 )“. 


+ ^l{l'+2)0^sClmsC(i-i)my + 2)(^ + 3) + (5;(J'_2)/3^sC(J+2)msC'(i+l)m 2 ~ 1)(^ ~ 2) + 0(y)(A7) 
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